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Frequency and phase synchronization of two coupled neurons with channel noise 
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We study the frequency and phase synchronization in two coupled identical and nonidentical 
neurons with channel noise. The occupation number method is used to model the neurons in the 
context of stochastic Hodgkin-Huxley model in which the strength of channel noise is represented 
by ion channel cluster size of neurons. It is shown that channel noise allows the two neurons to 
achieve both frequency and phase synchronization in the regime where the deterministic Hodgkin- 
Huxley neuron is unable to be excited. In particular, the identical channel noises lead to frequency 
synchronization in weak-coupling regime. However, if the coupling is strong, the two neurons could 
be frequency locked even though the channel noises are not identical. We also show that the relative 
phase of neurons displays profuse dynamical regimes under the combined action of coupling and 
channel noise. Those regimes are characterized by the distribution of the cyclic relative phase 
corresponding to antiphase locking, random switching between two or more states. Both qualitative 
and quantitative descriptions are applied to describe the transitions to perfect phase locking from 
no synchronization states. 
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I. INTRODUCTION 

The synchronization phenomena have been widely 
studied in neural systems in past decades [l|, 0, H, U, [^. 
Experiments show that the synchronization of coupled 
neurons could play a key role in the biological information 
communication of neural systems Q. Recent research 
also suggests that synchronization behavior is of great 
importance for signal encoding of ensembles of neurons. 
Especially the phase synchronization may be important 
in revealing communication pathways in brain [6] . Study- 
ing the synchronization of a pair of coupled neurons has 
attracted large amounts of research attention. In order 
to understand the dynamical properties of a neural net- 
work, it is important to characterize the relation between 
spike trains of two neurons in the network Q. What's 
more, studies show that noise enhances synchronization 
of neural oscillators. For example, the identical neurons 
which are not coupled or weakly coupled but subjected to 
a common noise may achieve complete synchronization. 
Actually, this is ageneral results for all the dynamical 
system d, i, [13, El, El- Both independent and corre- 
lated noises are found to enhance phase synchronization 
of two coupled chaotic oscillators below the synchroniza- 
tion threshold 

Among large population of neurons, different neurons 
are commonly connected to other group of neurons and 
receive signals from them. As a result of integration of 
many independent synaptic currents, those neurons re- 
ceive a common input sig nal which often approaches a 
Gaussian distribution 'J^. Therefore, noise was usually 
considered as external and introduced by adding to the 
input variables. However, recent work found that the 
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random ionic-current changes produced by probabilistic 
gating of ion channels, called channel noise or internal 
noise, also play an amazing role in single neuron's firing 
behavior and information processing progress [igI. ITtI. [l^ . 
Besides, Casado has showed that channel noise can allow 
the neurons to achieve both frequency and phase syn- 
chronization [20|, l2l| . This finding suggests that channel 
noise could play a role as promoter of synchronous neural 
activity in population of weakly coupled neurons. How- 
ever, Casado didn't give a quantitative description of the 
results. 

The magnitude of the ion channel noise is changed via 
the variation of the channel cluster size of neurons. It 
implies that synchronization in neural system is also re- 
stricted by the channel cluster size of neurons. Actually, 
the cluster size of ion channels embedded in the biomem- 
brane between the hillock and the first segment of neu- 
rons determines whether the neuron fires an action po- 
tential. The channel cluster size of this region is different 
for different neurons or for diff'erent developing stages of 
a neuron. With the decreasing of ion channel cluster size, 
ion channel noise would be induced thus the firing behav- 
ior of neurons would be greatly changed (for review see 
Refil7|). It is natural to ask to what extent the change 
of the channel cluster size affect the collective activities 
of neuron ensembles. For this purpose, we investigated 
the effect of ion channel cluster sizes (i.e, channel noise ) 
of neurons on synchronization of two coupled stochastic 
Hodgkin-Huxley (HH) neurons in this paper. 

Here we adopted a so called occupation number 
method rather than the Langevin method Casado had 
used to describe the single neuron for two reasons. First, 
the Langevin approach has been proved could not repro- 
duce accurate results for small and large cluster sizes. 
Second, occupation number method gives a direct rela- 
tion between channel cluster size of neuron and its firing 
behavior, and it's the fastest method for a given accuracy 
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[22l |. The main goal of our work is to explore what role 
the channel noise might play in the synchronization of 
two coupled neurons. We try to give qualitative as well 
as quantitative descriptions of the result. The practical 
meaning of our study is obvious, since the channel cluster 
size of neurons can be regulated by channel blocking ex- 
perimentally [23! , our study may provide a possible way 
to control neural synchronization. 

This paper is organized as follows. In section|ITl the oc- 
cupation number method of stochastic Hodgkin-Huxley 
neuron is introduced and the firing behaviors of neurons 
with different channel cluster sizes are demonstrated. In 
the following sections, we explored the combined effect 
of coupling strength and cluster size on the synchroniza- 
tion behaviors of two neurons with an electrical synap- 
tic connection. Section IIIII is devoted to frequency syn- 
chronization. The phase synchronization of identical and 
nonidentical neurons are discussed in Section ITVl A con- 
clusion is presented in Section IVl 



II. THE MODEL 



TABLE I: Parameters and Rate Functions Used in Simula- 
tions. 
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Hodgkin-Huxley neuron model provides direct relation 
between the microscopic properties of ion channel and 
the macroscopic behavior of nerve membrane. The mem- 
brane dynamics of HH equations is given by 
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where V is the membrane potential 
Vl are the reversal potentials of the potassium, sodium, 
and leakage currents, respectively. Gk, G^a, and Gl are 
the corresponding specific ion conductances. C™ is the 
specific membrane capacitance, and / is the current in- 
jected into this membrane patch. The voltage-dependent 
conductances for the and Na^ channels are given by 
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where N'^'^^ 



and iVX" 



are the numbers of open potas- 
sium and sodium channels. S is the membrane patch 
area, and ^Na give the single-channel conductances 
of if + and 7Va+ channels. Then the numbers of total 
potassium and sodium channels Nk and iV^va are given 
by the equations Nk = Pk x S and N^a = PNa^S, where 
Pk and p^a are the and Na^ channel densities re- 
spectively. By introducing time constants tk = ^^"^ , 

TNn = — — and Tr = , we end up with the foUow- 
ing equation for the membrane potential 
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Individual channels open and close randomly. If the 
number of channels are large and they act independently 
of each other, then, from the law of large numbers, 
jyopen^ iVif ( or N'^^^ /N^a) is approximately equal to 
the probability that any one K'^ (or iVa^) channel is in 
an open state, and can be represented as continuous de- 
terministic gating variables and m^h. This leads to 
the deterministic version of HH model 124, Uffl , 
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where — px x and 'gj^a — PNa x jNa are the max- 
imal potassium and sodium conductance per unit area, 
n"* indicates that the channel has four separate gates 
and that a channel is opened when all those gates 
are open; m^h indicates that three m-gates and one h- 
gate must be opened to open a A^a+ channel. The gating 
variables obey the following equations. 



x) — l3x{V)x,x = m,h,n, (5) 



where a^iV) and l3x{V) {x — m, ft,, n) are voltage depen- 
dent opening and closing rates and are given in Table |l] 
with other parameters used in the simulations. 

The deterministic HH neuron model [Eqs. ([4])- ([5])] de- 
scribes the transmembrane potential without the need 
to treat the underlie activity of individual ion channels. 
However, for the limited number of channels, Eq. @ is 
no longer valid and statistical fluctuations will play a role 
in neuronal dynamics [l^. So we have to return to Eq. 
^ and have to determine N'^'^^ and N"^^^ as a function 
of time by stochastic simulation methods based on state 
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FIG. 1: Kinetic scheme for a stochastic potassium channel (a) 
and sodium channel (b). n4 and m^hi are open states, while 
the other states are no-conducting. 



diagrams that indicate the possible conformation states 
of channel molecules. 

As shown in Fig. [1] both and Na^ channels exist 
in many different states and switch between them accord- 
ing to voltage depended transition rates (identical to the 
original HH rate functions), [n.i] is the state of K'^ chan- 
nel with i open gates and [rriihj] is the state of iVa+ 
channel with i open m-gates and j open /i-gates. Hence, 
[n4] labels the single open state of the channel and 
[ma/ii] the Na^ channel. Usually we can simulate the 
kinetic scheme of each ion channel to get the numbers 
of open sodium and potassium channels at each instant. 
However, it is not an efficient way because many tran- 
sitions of states do not change the conductance of the 
channel. Instead of keeping track of the state of each 
channel, we keep track of the total populations of chan- 
nels in each possible state so N'^'^"' and N'^^"' at each 
instant can simply be determined by counting the num- 
bers of channels in state [n^] and [ma/ii]. Specifically, 
if the transition rate between state A and state B be 
r and the number of channels in these states be ua and 
ub- Then, the probability that a channel switches within 
the time interval {t,t + At) from state A to S is given by 
p = rAt. Hence, for each time step, we determine Auab, 
the number of channels switch from A to B, by choosing 
a random number from a binomial distribution [iSl. 
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Then we update ua with ha — Ahab, and hb with ub + 
AuAB ■ To make sure that the number of channels in each 
state is positive, we update those number sequentially, 
starting with the process with the largest rate and so 
forth. 

Voltage-gated ion channels are stochastic devices. The 
origin of channel noise is basically due to fluctuations 
of the fraction of open ion channels (thus the channel 



FIG. 2: The mean firing frequency as a function of the channel 
cluster size A?^ for 7 = nA/cm^{o), 1 = 3 iJ,A/cm'^{n) and 
1 = 6 fxA/cm?{A). The data are obtained from spike trains 
of 2000 action potentials. 



currents) around the corresponding mean values. The 
strength of the fluctuation is inversely proportional to 
the number of total ion channels [13, [H. Though aver- 
age membrane current is at constant, as we will see, the 
variance of the Na~^ and currents cause remarkable 
effects on neuronal dynamics. In this work, we introduce 
channel cluster size N {N = Nk = A^Afa/S) as a mea- 
surement of channel noise level so that the correct propor- 
tion between Na'^ and channel densities is preserved. 
With increasing channel cluster size, the fluctuations of 
the fraction of open ion channels, thus the variance of 
the the corresponding channel currents decreases. For 
a large number of channels this noise becomes negligible 
(i.e., the deterministic case). The threshold constant cur- 
rent for deterministic HH neuron to generate consecutive 
action potentials is Ith = 6.26fiA/crn^. However, due to 
the channel noise, the stochastic HH neurons can gener- 
ate spiking activity with subthreshold input current . 
Fig. [2] shows the mean firing frequency (defined in Sec. 
imp as a function of channel cluster size for different con- 
stant current. If the channel cluster size is small, the 
neuron fires action potentials with high frequency. As 
the channel cluster size increases, the mean firing fre- 
quency drops quickly, approaching the deterministic case 
that no firing activities occur with the same subthreshold 
input currents. With decreasing the input current, the 
firing frequency decreases. However, the firing activities 
will not vanish if the input current is decreased to zero. 
Thus, as have demonstrated, channel noise shifts the on- 
set of firing behavior to lower values of input current /. 

To explore the synchronization phenomena, we con- 
sider two stochastic HH neurons coupled by an electrical 
synaptic connection. The system is described by the fol- 
lowing equations. 
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FIG. 3: The winding number u\/uo2 as a function of N1/N2 
with iVa = 2 X 10^ and e = 0.3 (o), e = 0.1 (•), e = 0.01 (*), 
£ = 0.001 (□). 



Here Vi and V2 are the instantaneous membrane poten- 
tials of the two neurons and e is the diffusive coupHng 
strength between the neurons. N'^^^, -/V^^^", N'^^"' , 
^atqT numbers of open and Na^ channels of 

neuron 1 and 2, respectively ; Nk^ N^an Nk2, -^Wag are 
the numbers of total and iVa+ channels for neuron 
1 and 2, respectively. Ii and I2 are two constant input 
currents which are set at h = I2 = QjiA/crrP'. Here, A^i 
and N2 [Ni — — N^ai/i, « = 1,2) are the channel 
cluster sizes for each neuron. 

The numerical integration of system mentioned above 
is carried out by using the Euler algorithm with a step 
size of 0.01ms. And all simulations are working in Ito 
framework. The occurrences of action potentials are de- 
termined by upward crossings of the membrane potential 
at a certain detection threshold of lOmV if it has previ- 
ously crossed the reset value of -50mV from below. 



III. FREQUENCY SYNCHRONIZATION 

From the above-mentioned system, we can get two 
point processes in the following form. 



N 



(9) 



Each one gives the spike sequence of a particular neuron. 
The mean spiking frequency of neuron i (i = 1,2) is 
defined as. 
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Generally, synchronization means an adjustment of 
timescales of oscillations in systems due to their cir- 
cumstances. In other words, oscillators can shift the 
timescales to make their ratio close to a rational number 
n : m, where n and m are integers. This phenomenon 
is usually referred to as n : m frequency synchroniza- 
tion, and its suitable measure is the closeness of the ratio 
of u>i/uj2 to the chosen rational number n : m [25]. In 
this paper we will discuss only 1 : 1 synchronization. 
Note that the frequency locking discussed in this section 
is in a stochastic sense, and refers to the equivalence of 
the average frequencies rather than the instantaneous fre- 
quency. So it is not a sufficient condition for synchroniza- 
tion. However, since the firing rate of a neuron is often 
argued to carry information of the stimulus, studying of 
the frequency synchronization is especially meaningful in 
the context of rate coding scheme. 

We investigated the shift of winding number loi/uj2 
along with both the variation of coupling strength and 
channel cluster size. When the coupling strength is 
small, as shown in Fig. [31 with the increasing of N1/N2 
from N1/N2 < 1, the winding number will decrease 
from a value at which the two neurons are not synchro- 
nized. When two channel cluster sizes are the same ( 
N1/N2 = 1), both neuron will be frequency locked. Fur- 
ther increasing of N1/N2 will desynchronize them to a 
certain levels. Note that in the range of N1/N2 > 1, in- 
creasing the value of coupling strength tends to increase 
LUi/uj2, and the increasing is larger with larger coupling 
strength. In panel A of Fig. 21 it is seen clearly that 
though frequency synchronization can be achieved with 
arbitrary chosen value of A*'2, the tuning is very critical, 
as a small variation of Ni from A^2 leads to desynchro- 
nization. As has been described above, for a isolated 
neuron, its average firing rate decreases with increasing 
its channel cluster size (i.e., decreasing the channel noise 
intensity). Thus, at a given channel cluster size, the chan- 
nel noise term is identical for the two neurons. In this 
case, their average firing rates would be the same, which 
means the two neurons are frequency locked. 

However, if the coupling strength is increased to a 
rather large value (for example, e = 0.3 in Fig. [3]), the 
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FIG. 4: The winding number iOi/uj2 as a function of N1/N2 with e = 0.01 (A), e = 0.1 (B) and different iV2: iV2 = 2 x 10^ (o), 
iV2 = 2 X 10^ (□) and iVa = 2 x 10* (*). 



coupling strength starts to take command of the fre- 
quency synchronization as the two neurons are able to 
be entrained in a wider range of channel cluster size (i.e., 
channel noise level). It is seen in panel B of Fig. |4] that 
neurons with larger value of N2 is easier to get frequency 
entrained in a wider range with lower coupling strength. 
Even in the weak-coupling shown in panel A of 

Fig. m large channel cluster sizes tend to enhance syn- 
chronization (make uji/uj2 closer to 1). This implies that 
neurons with large channel cluster sizes (i.e., small chan- 
nel noise level) are easier to adjust their timescales to 
make their firing rates close to each other. However, if 
the channel noises are too small, the neurons won't fire 
spikes with subthreshold stimuli. 

It is concluded that for identical, symmetrically cou- 
pled neurons, when the coupling strength is small, the 
channel cluster sizes at frequency synchronization must 
be the same, whereas the coupling strength only has 
a limited effect only when the channel cluster sizes of 
the two neurons are not same. However, when the cou- 
pling strength is rather large, the two neurons are able to 
achieve frequency synchronization with a greater range of 
channel cluster sizes. In this regime, though large chan- 
nel noise degrade frequency synchronization, small chan- 
nel noise intensities help to get the neurons frequency 
synchronized with subthreshold stimuli. 
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FIG. 5: Synchronization diagram for the distribution of 
cyclic relative phase P{^)- Region 1 correspond to state of 
monomodal distribution (•), region 2 to bimodal distribution 
(o), region 3 to drifting evolution of Ai^ (□); region 4 to no 
firing area (x). There are no lines plotted to separate the 
subregions in those regions because the distribution of -?($) 
changes in a continuous manner in those regions. 



IV. PHASE SYNCHRONIZATION 



cations like neural spike sequence, it is useful to define 
an instantaneous phase (f)(t) by linear interpolation. 



Given a data set or some model dynamics there exists a 
variety of methods to define an instantaneous phase (j){t) 
of a signal or a dynamics [131 . However, for a stochastic 
system it is essential to assess the robustness of the phase 
definition with respect to noise. In many practical appli- 



0(0 = 27r ^ ^" + 27m (i„ ^ t sC t„+i), (11) 

where i„ is the time at which the neuron fires a spike. 
The instantaneous phase difference between them is then 
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FIG. 6: The distribution of the cychc relative phase -P($) corresponding to some representative points of the synchronization 
diagram in Fig. \5\ (A) e = 0.3, iV = 2 x 10^ (B) e = 0.3, TV = 2 x 10''; (C) e = 0.3, N= 2 x 10^ (D) e = 0.08, = 2 x 10^ 
(E) e = 0.08, A = 2 x 10^ (F) e = 0.08, A = 2 x 10^; (G) e = 0.02, A = 2 x 10^ (H) e = 0.02, A = 2 x 10^ (I) e = 0.004, 
A = 2 x 10^. Each plane have different vertical scales. 



given by 



(12) 



Phase synchronization is a weak form of synchronization 
in which there is a bounded phase difference of two sig- 
nals. Usually, the relative phase can vary from — oo to 
+00 in stochastic system if the coupling is weak and/or 
the noise level is high. However, if we increase coupling 
strength and adjust noise to a low level, the relative phase 
will fluctuate around some constant values. Sometimes, 
noise would induce a phase slip where the relative phase 
changes abruptly by ±27r. Thus, it is useful to define the 
phase locking condition in a statistical sense by the cyclic 
relative phase T| 



$ = tp{mod2n) 



(13) 



A dominant peak of the distribution of this cyclic relative 
phase P{^) reflects the existence of a preferred relative 
phase for the firing of both neurons. When this pre- 
ferred phase is zero we speak of phase synchronization in 



a statistic sense. We speak of out-of-phase synchroniza- 
tion when the distribution P($) peaks around a nonzero 
value of Especially if the nonzero value is tt, we call 
it antiphase synchronization [U . 



A. phase synchronization of identical neurons 

In this section, we study phase synchronization of two 
identical neurons (A^i — N2 ^ N). In Fig. we present 
the synchronization diagram in terms of the coupling 
strength e and channel cluster size N. A different form 
of the distribution P{^), which is plotted in Fig.[6l char- 
acterizes each region in it. The corresponding temporal 
evolution of the relative phase is illustrated in Fig. [T] 
We will give a detailed description of each region in the 
following part. 

In region 1, the distribution shows a monomodal char- 
acter as plotted in panel A, B, C of Fig. [S] In this re- 
gion, when both e and N are very large, the distribu- 
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FIG. 7: Temporal evolution of the relative phase (A) e — 0.3, 
N = 2x10^ (B) e = 0.3, iV = 2x10^ (C) e = 0.08, N = lO''; 
(D) e = 0.02, iV = 2 X 10^ 

tion of P{^) is a narrow peak on tt. Thus we can speak 
of antiphase synchronization for a statistic out-of-phase 
locking is achieved. With the decreasing of N, the peak 
is still on TT but becomes broader(see the change of peaks 
in A ^ i3 — » C ). This suggests that in the case of 
large coupHng strength, large channel cluster size (i.e., 
small channel noise) allow the statistical antiphase syn- 
chronization to approach full antiphase synchronization 
appearing in deterministic systems. As can be seen in 
Fig. [5l there is a minimal value of the coupling strength 
£o for which the antiphase locking becomes stable in a 
statistical sense. This minimal coupling strength Eq (~ 
0.115) is independent of the channel cluster size. To show 
the effect of channel noise on antiphase synchronization, 
we demonstrated the temporal evolution of relative phase 
in this situation in panel A and B of Fig. [T] Obviously, 
decreasing channel cluster sizes will lead to larger fluc- 
tuation of the relative phase due to channel noise, thus 
give distribution of P($) a broader peak, but it does not 
destruct antiphase synchronization. 

Region 2 marked by open circles corresponds to the 
bimodal distribution of P{^), as shown in panel D, E 
and F of Fig. [51 In this region, when e and N are large 
[2(a)], the two peaks of the distribution are well spaced 
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FIG. 8: Synchronization indices for two identical neurons ver- 
sus coupling strength with various ion channel cluster sizes A''. 



from each other. To uncover the underlying mechanism 
of this bimodal distribution, we investigated the tempo- 
ral evolution of relative phase in this situation. As we 
observed in panel C of Fig. [71 '0(0 "^ill fluctuates suc- 
cessively around one of a pair of symmetric values for a 
long period then switch suddenly to the other one. This 
fact clearly reflects the two-state character of the phase 
dynamics. We argue that this two-state dynamics is the 
result of a compromise between coupling and noise. The 
existence of a two-state dynamics suggests the possibil- 
ity of inducing a kind of stochastic resonant behavior 
by coupling the relative phase to the action of an ex- 
ternal periodic forcing [21]. Again, if we decrease N to 
enter 2(b) area, the two peaks will be broader, and grad- 
ually overlapped and have small wings. The overlapping 
means that the bimodal distribution becomes unstable. 
The wings implies that preferred relative phase for the 
firing of both neurons does not prefer to some certain 
values anymore, thus the synchronization becomes weak. 
If we further decrease N to enter 2(c) area, we find that 
the two peaks move closer and then merge but still have 
a maxima around tt (see panel F of Fig. [6l) 

Actually, the phase-locking phenomena in noise-free 
neural systems have been well studied through effec- 
tive coupling analysis [H, [H, [s^]- S. K. Han and Ku- 
ramoto had demonstrated that diffusive interaction will 
dephase interacting oscillators and may stabilize them 
at a phase difference given by the corresponding stable 
fixed point according to the initial condition (see detail 
in Ref. J2§\). In our case, those stable fixed points are 
stochastic variables with single peaks distribution alike 
the peaks demonstrated in Fig. [SI In region 1, there 
is only one stable fixed point distributed around tt and 
will become broader if the noise intensity is increased. 
In region 2, the system has two stable fixed points. The 
system will be stabilized at one point according to the ini- 
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tial condition, then the channel noise occasionally change 
the initial condition and stabilized the system at the 
other one. If the channel cluster size N is extremely 
large, though there are still two stable fixed point, the 
channel noise is too weak to change the initial condition 
frequently, and only one of the two peaks can be ob- 
served with certain recording time interval(not shown). 
As N decreases, the channel noise becomes larger, giving 
broader distributions of the two stable fixed points and 
more frequent switches between them. It is the broader 
distributions of the two stable fixed points that leads the 
overlapping of the two peaks in this area. 

In region 3, we find the bimodal distribution of -P(<&) 
mentioned before disappears and the single peak distri- 
bution appears again with large wings. Panel G of Fig. [6] 
shows a representative cyclic relative phase distribution 
corresponding to the 3(a) area, which is characterized by 
a peak around n and another smaller one around 0. If we 
decreasing the coupling strength to enter 3(b) area from 
2(c) area, the central peaks will decrease in height [panel 
H ] and eventually disappear [3(c) area, panel I]. In 3(c) 
area, the relative phase will drift unboundedly and P(<&) 
ceases to be useful (see panel D of Fig. [T]). When the 
coupling strengths is very weak, the system at hand can 
be considered as two uncoupled neurons which fire inde- 
pendently due to ion channel noise. Thus their relative 
phase can be at any arbitrary value (as show in panel D 
of Fig. [7]), and gives relative phase a smooth distribution. 

Region 4 in Fig. [5] is the silent state in which both 
neurons cannot fire spikes but only perform low ampli- 
tude fluctuations around its resting potentials under the 
combining effects of coupling and channel noise. 

Next, we characterize those peaks with synchroniza- 
tion indices which are defined as 



7^ = (cos$(t))2 -I- (sm$(t))^ 



(14) 



where (...) denotes temporal averaging. The index 7 as- 
sumes values between (no synchronization) and 1 (per- 
fect phase locking) [25i] . 

Fig.[5]quantitatively demonstrated the synchronization 
of two identical neurons under the effect of both coupling 
strength and channel noise. When coupling strength is 
small, the two neurons show almost no synchronization 
[7 w 0, corresponding to 3(c) area of Fig. [5] or silent 
state when N is large(incomplete lines, corresponding 
to region 4 of Fig. [5]). The synchronization index 7 is 
not sensitive to the change of channel cluster size N in 
small coupling strength region. As the coupling strength 
increases the degree of synchrony of the two neurons in- 
creases. The maximal synchrony appears at 7 « 1 (corre- 
sponding to region 1 of Fig. [5]). Note that when channel 
cluster sizes are large, as shown in Fig. [SJ there exists 
a step-like transition (a threshold) to perfect synchro- 
nization. The threshold is around about e = 0.04 when 
N = 2 X 10^. However, it disappears when cluster sizes 
decrease, and the transition becomes a graded type. It 
implies that channel noise can 'soft' the threshold to give 
a wider range of synchronization degree. 




1 2 A 6 2 4 6 



CD (Rad) 



FIG. 9: The distribution of the cyclic relative phase P{^) for 
nonidentical neurons. (A)e = 0.2, iVi = 2 x 10^, N2 =2x 10^; 
(B)e = 0.08, iVi = 1 X 10^ iV2 = 2 x 10''; (C)e = 0.03, A^i = 
2 X 10^, = 2 X 10^. Each plane have different vertical 
scales. 
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FIG. 10: Synchronization indices for two nonidentical neurons 
versus coupling strength with channel cluster size of the first 
neuron = 2 x 10^. 



B. phase synchronization of nonidentical neurons 

Actually, neurons in nature are not identical. The non- 
identity can be achieved in numerical simulation by mis- 
matching neuronal parameters (the leakage conductance 
gi, for example) |15j . Here, we introduce a mismatch 
into the channel cluster sizes of the two neurons (i.e., 
Ni ^ iV2). With this parameter heterogeneity, the neu- 
ron with smaller channel cluster size, due to its larger 
channel noise, is easier to be excited by subthreshold 
stimuli and has larger firing rate than another one. 

In the case of two nonidentical neurons, the above 
mentioned cyclic relative phase distributions are still ten- 
able and the perfect phase synchronization can also be 
achieved (see Fig. [TU]) . However, there are three excep- 
tions. First, because the symmetry of the distribution 
P($) is dependent on the symmetry of the system, for 
nonidentical neurons the distribution of the cyclic rel- 
ative phase is asymmetric (see Fig. [5]). This fact was 



9 



also confirmed by applying different tonic subthreshold 
currents to two neurons to make the system asymmetric 
pll |. Second, as mentioned before, the two weakly cou- 
pled identical neurons with large channel cluster sizes are 
unable to fire spikes under a subthreshold stimulus (see 
the plot of = 2 X 10^ in Fig. [8|). However, when a neu- 
ron with large channel cluster size coupled with a small 
one which could be excited by subthreshold stimulus due 
to channel noise, the large one is excited by the small 
one through coupling. As shown in Fig. (TUl comparing 
with identical situation, the neuron can be excited when 
N — 2 X 10^ and e is rather small. It is also seen that in 
the weak coupling region, the identical neurons exhibit 
higher degree of synchronization than the nonidentical 
ones. Whereas in the strong coupling region, when a neu- 
ron is coupled to another one which has larger N, they 
exhibit higher degree of synchronization. This is consis- 
tent with the frequency synchronization case where iden- 
tical neuron is frequency locked even coupling is weak, 
but nonidentical neuron can also be synchronized if the 
coupling is strong, and neurons with large channel clus- 
ters are easier to get synchronized. 

V. CONCLUSION 

In conclusion, the frequency and phase synchronization 
of two coupled stochastic Hodgkin-Huxley neurons are 
studied by varying coupling strength and channel clus- 
ter sizes. The two neuron is coupled via a gap junction 
because the gap-junctional (diffusive) coupling can gen- 
erate rich dynamical behavior [3ll|. What's more, with 
this simple coupling, we could emphasized on the effects 
of channel noise and ignore the inessential details of com- 
plex synaptic process. Our studies show that when the 
coupling is weak, the cluster sizes of the two neurons 
must be the same to achieve frequency synchronization, 
and the synchronization region is very narrow. How- 



ever, when coupling is strong, the two neurons can be 
frequency entrained in a wide region. For two identical 
neurons, a state of statistical antiphase synchronization 
is reached in the strong coupling region if the cluster 
size is large enough. In this state, the relative phase be- 
tween two spike trains would be around tt. As the cou- 
pling strength and channel cluster size are reduced, the 
phase-locking condition is lost and a rather complex be- 
havior would appear. This complex behavior, as we have 
argued, is the result of a compromise between coupling 
strength and channel noise. We use synchronization in- 
dices to characterize the transitions to synchronization, 
and find that there exit a threshold to synchronization. 
When channel cluster sizes are small, channel noise can 
'soft' this threshold to present synchronization at a wider 
range. For two nonidentical neurons, the distribution 
of the cyclic relative phase is asymmetric and the silent 
state in identical situation disappears. This, as we have 
pointed out, is due to asymmetric of the system and spon- 
taneous firing induced by channel noise. 

It is helpful that our study is important for the un- 
derstanding of coupled stochastic systems and possible 
applications especially in neuroscience where the syn- 
chronization activity could be tuned through the control 
of the channel noise via channel blocking. By applying 
channel cluster size control to real neural systems one 
should be able to influence neural synchrony. Further 
work should focus on more sophisticated models and on 
coupling more than two neurons. 
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